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Abstract. Action of solar wind on arbitrarily shaped interplanetary dust particle is investigated. 
The final relativistically covariant equation of motion of the particle contains both orbital evo- 
lution and change of particle's mass. Non-radial solar wind velocity vector is also included. The 
covariant equation of motion reduces to the Poynting-Robertson effect in the limiting case when 
spherical particle is treated, the speed of the incident solar wind corpuscles tends to the speed 
of light and the corpuscles spread radially from the Sun. The results of quantum mechanics have 
to be incorporated into the physical considerations, in order to obtain the limiting case. The 
condition for the solar wind effect on motion of spherical interplanetary dust particle is p'out ~ 
( 1 — o-'pr/o'tot ) Pini where Pi„ and p'out a-r^ incoming and outgoing radiation momenta (per 
unit time) measured in the proper frame of reference of the particle; a'p^ and a't^t are solar wind 
pressure and total scattering cross sections. 

Analytical solution of the derived equation of motion yields qualitative behaviour consistent 
with numerical calculation. This holds also if decrease of particle's mass is important, when also 
outspiralling from the Sun may occur. 

Real flux density of radial solar wind energy yields that the time of spiralling toward the Sun 
may differ in about 10 % from the case of constant flux for heliocentric distances smaller than 
about 10 AU (if radius of the particle is unchaged). The differences between the real solar wind 
effect (when also non-radial solar wind velocity vector is taken into account) and the standard 
approach may be even more significant for heliocentric distances greater than 10 AU: solar wind 
can cause outspiralling from the Sun, even if radius of the particle is constant. 

Real fiux density of solar wind energy produces shift of perihelion of interplanetary dust 
particles. This result significantly differs from the standard treatment of the action of the solar 
wind on dust particles, when analogy with the Poynting-Robertson effect is stressed. Moreover, 
the evolution of the shift of perihelion depends on orbital position of the parent body at the 
time of ejection of the particle. 

Key words, cosmic dust, electromagnetic radiation, solar wind, relativity theory, quantum theory, 
equation of motion, orbital evolution 

1. Introduction 

The Poynting-Robertson effect is used in modelling of orbital evolution of dust grains under the action 
of electromagnetic radiation (of the central star), for many decades (e. g., Foynting 1903, Robertson 
1937, Wyatt and Whipple 1950, Dohnanyi 1978, Kapisinsky 1984, Jackson and Zook 1989, Leinert 
and Griin 1990, Gustafson 1994, Dermott et al. 1994, Reach et al. 1995). It was presented that solar 
wind operates in a similar way. The action of the solar wind on motion of interplanetary dust particle 
was discussed, in a heuristic way, e. g., by Whipple (1955). He has mentioned also the results of 
laboratory experiments: intense bombardment of a material by energetic corpuscles destructs the 
material and this effect is known as a " sputtering" . Current opinion is that we have two different 
effects of the solar wind: i) motion of dust particle is influenced by the incident solar wind, and, ii) the 
corpuscular sputtering decreases mass of the particle (e. g., Whipple 1955, Dohnanyi 1978, Kapisinsky 
1984, Leinert and Griin 1990). There has been an attempt to better understand physics of the action 
of the solar wind on the motion of dust particle. As a first attempt, we can mention Robertson and 
Noonan (1968, pp. 122-123). The authors formulated relativistically covariant equation of motion of 
the particle under the action of the solar wind. However, their result does not admit any destruction of 
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the particle. A more realistic view was presented in Klacka and Saniga (1993), where also space-time 
formulation of the problem was suggested. As a result, corpusciilar sputtering is an indispensable part 
of the equation of motion for the action of solar wind on interplanetary dust particle. 

Our paper presents space-time formulation of the action of solar wind on arbitrarily shaped inter- 
planetary dust particle. Equation of motion in a relativistically covariant form is derived. Moreover, in 
order to be physically correct, results of quantum theory are also taken into account. The results of the 
paper are consistent with the results of the papers by Klacka (2008a, 2008b) for the electromagnetic 
radiation. Our theoretical derivations hold for any solar wind velocity vector and the result can be 
easily applied to other stars with stellar winds. 

Application of the derived equation of motion to spherical interplanetary dust particle is presented 
in the form of orbital evolution of the particle. While submicron dust particles are driven mainly 
by Lorentz force (motion of charged particles in the interplanetary magnetic field, Dohnanyi 1978, 
Leinert and Griin 1990, Dermott et al. 2001), collisions among particles are important for particles 
of radii larger than hundred of micrometres, approximately (Griin et al. 1985, Dermott et al. 2001). 
This paper deals with orbital evolution of micron-sized spherical interplanetary dust particles, when 
effects of solar gravity, solar electromagnetic radiation and solar wind (solar corpuscular radiation) 
are relevant. Radial solar wind is conventionally used. However, the newest observations ( Bruno et 
al. 2003) show that velocity vector of the solar wind corpuscles is non-radial and the angle between 
the velocity vector and radial direction is practically independent on heliocentric distance. Our paper 
compares orbital evolution of spherical interplanetary dust particle for standardly used approach of 
time independent radial solar wind and more real solar wind model. Also the action of mass decrease 
of the interplanetary dust particle is taken into account. 

Section 2 derives relativistically covariant equation of motion of an arbitrarily shaped particle under 
the action of solar wind (including non-radial component of the solar wind velocity). Sec. 3 summarizes 
important equations for the Poynting-Robertson effect. Equation of motion of the spherical particle 
under the action of solar corpuscular and electromagnetic radiation, and, solar gravity, is given in Sec. 
4. Sec. 5 deals with secular evolution of particle's orbital elements under the action of solar radiation 
(electromagnetic and corpuscular - solar wind), in an analytical way. The following section Sec. 6 
concentrates on detail treatment of the numerical results and compares the results for conventional 
time independent radial solar wind with those obtained for more real solar wind model. 



2. Equation of motion - solar wind effect 

Relativistically covariant equation of motion of interplanetary dust particle is derived in this section. 
Our derivation enables understand physics of the action of the solar wind on the motion of dust 
particle (compare heuristic explanation by Whipple 1955 and space-time formulations presented by 
Robertson and Noonan 1968, pp. 122-123, and, Klacka and Saniga 1993). 

We will show that the corpuscular sputtering is an indispensable part of the action of the solar 
wind on the interplanetary dust particle, so the sputtering cannot be considered as an another effect 
of the solar wind. Moreover, non-radial solar wind velocity vector can be easily incorporated into the 
final equation of motion. Finally, the covariant formulation yields the Poynting-Robertson effect in 
the limiting case when speed of the solar wind corpuscles tends to the speed of light. The limiting 
case is fulfilled under the assumption that the total cross section of the interaction between the solar 
wind corpuscles and the interplanetary dust particle is given by the results of quantum theory and 
not by the classical non-quantum physics. 

This section presents also results with an accuracy to the order (v/u)^, where v is orbital velocity 
of the particle with respect to the Sun and u is the solar wind speed. These results will be used in 
practical modelling of orbital evolution of interplanetary dust grains in Sees. 4 and 5. 
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2.1. Incident radiation 

Let us introduce two inertial reference frames. The first is the proper reference frame of a particle 
moving with velocity v around the Sun. The particle is at rest in it's proper frame of refc^rence. 
Quantities measured in this frame will be primed. The second frame is associated with the Sun. This 
frame is "stationary reference frame". 

We will suppose that all corpuscles of the solar wind are of the same mass mi and of the same 
velocity u (or u' in the proper reference frame of the particle). Thus, each of the corpuscles has the 
following four-momentum 

p'^,={E[/c; p\) = m,^(u')(c; u') (1) 

in the proper reference frame of the interplanetary dust particle, or 

= {E^/c ; pi) = mi7 (u) (c ; u) (2) 

in the stationary reference frame; c is the speed of light. 

Let a beam of such solar wind corpuscles hits the dust grain. Energy and momentum incident on 
the particle per unit time in its proper frame are 

E'in = ^tot n' u' E[ , 

P'in = c^'tot n' v! p'l , (3) 

where n' is the concentration of the solar wind corpuscles and a'^g^ is the total scattering cross section 
of the interplanetary dust particle. Using Eq. (1), we can rewrite Eqs. (3) into the form of the incident 
four-momentum per unit time 

( E' 

P'tn = '^tot n' u' ( ; p'l 



l,;..„'„'i5;(i;^) . 



(4) 



Introducing the flux density of the incident energy (energy flow per unit area perpendicular to the 
beam of solar wind corpuscles per unit time) 

S' = n'u'E[ , (5) 
Eq. (4) can be rewritten to the form 

/M _ (Ek ■ w \ 

pTn = I s' al, (i ; ^) • (6) 

Having a four-vector B''^ = {B'^ ; B') in the proper reference frame, the components of the 
four-vector in the stationary reference frame are given by generalized special Lorentz transformation: 



S° =7(^)(b'° + ^) , 
B = B' + ([7(.)-l]^ + ^-^bAv 



(7) 



or inverse 
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Now, using Eqs. (6) and (7) we get 



(8) 



Pi', 



Pin 



~2 ^ ^tot 



U 



V ■ u 



(9) 



We have to express the primed quantities (except of (r'tot) the right-hand sides of Eqs. (9), i.e. 5' 
= n'u' E'l and u' , through the unprimed quantities measured in the stationary reference frame of the 
Sun. The energy E'l we obtain from Lorentz transformation of to the proper reference frame of the 
interplanetary dust particle. It holds 

1 E^= uj E^, (10) 

where we defined the quantity 



(11) 



Other quantities we get from transformation of the four- vector of the current density = (nc ; nu) 
to the corresponding four- vector j"^ = [n'c ; n'u'). The transformation yields 



n = u n , 

, 1 

u = — a , 

CO 



where the vector 
a = u + 
has magnitude 



a = < u + 7 (v) V - 27 (u) t) • M + 7 (v) 



V ■ u 



1/2 



Thus, u' = a/io and the flux density of energy is 
S' = ^S, 



s 



u 

nuEi 



(12) 
(13) 

(14) 
(15) 



according to Eqs. (5), (10) and (12). 

Finally, using Eqs. (9), (11), (12), (13) and (15), one obtains 



Pin ~ ^tot * 



Pi 



a CO 1 u 



a. 



tot 



C U CO c 

The incident four-momentum of solar wind per unit time is 

1 . _ a CO .. 



(16) 
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2.2. The reaction of spherical dust particle on the incident solar wind 

The incident solar wind corpuscule may be reflected from the surface of the interplanetary dust particle 
(IDP) or may cause its erosion/destruction and decrease the mass of the IDP, or. in general, similar 
processes as reflection, absorption and diffraction may occur. Let the particle's loss of energy (per unit 
time) in the proper reference frame of the IDP, be E'^^^ — E'^^ . iJ^^j can be written as an a;'-part of 
the incident energy per unit time. The relation 



E' = r' F/. 



(18) 



holds for the outgoing energy. In order to express the outgoing monicintuin wc; dciclarc; the ortonormal 
vector basis {f'j]j = 1,2,3} in the proper reference frame of the particle and three velocity vectors 
{u'j = u' f'j; j = 1,2,3} corresponding to these unit vectors. We suppose that u'l = u' . Now, the 
outgoing momentum per unit time is 



P'out 



1 



(19) 



where (jp,. ^ (j = 1, 2, 3; tTp^ ^ = a'^^ are pressure cross sections (analogy with optics - see Klacka 
2008a, 2068b). 

The outgoing four-momentum per unit time is given by Eqs. (6), and (18)-(19): 



I' out — 



j P out 



P out 



'tot 



S' x' 



<t) "^'^ c c f • 



(20) 



Generalized special Lorentz transformation of using Eqs. (15), gives for the outgoing four- 

momentum per unit time in the stationary reference frame 



M _ 1 / a w , 

Pout — ~ '^tot '-' ^ 

C U C 



+ 



tot 



1 



^tot S 



''tot - — 2^ ^ 



C 

3 _/ 



j=2 ^tot 



a u! 
u 



TJP. 

■' c 



c 



where 

= {j{v) c ; 7(w) v) 
is four- velocity of the IDP. The other four- vectors are 

1 1 Ui 



= 

^3 

, , / V ■ u'j 



u 



u 



{^{v)-l) 



V ■ u 



J = 1,2,3, 
and oji = = Ml = u. 



(21) 



(22) 



(23) 
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2.3. Equation of motion 



Now we can write equation of motion of the IDP under the action of the solar wind, in a relativistically 
covariant form: 



dT 



Pout 



(24) 



Eq. (24) yields, using Eqs. (17) and (21), 



d^ 



= - a. 



S 



a (jO 



H 




x' - 












V ^totJ. 





3 / 

pr,j 



j=2 



(25) 



where = m U^^ is four-momentum of the IDP of mass m and r is the proper time of the particle. 
Using 



dp^' d , dm dt/^ 

ar ar ar ar 



(26) 



Eq. (25) yields not only acceleration of the particle, but also change of the particle's (rest) mass, due 
to the interaction of the IDP with the solar wind. The change of the mass is given, on the basis of 
Eqs. (11), (17), (22), (25) and (26), by the expression {U^JJ^ = c^, f/^ dU^'/dT = 0) 



— - - - a' S—(x'-l) 



(27) 



One can easily verify that Eq. (27) corresponds to the famous Einstein's equation dm/dT = 

iE'i„ - E'^^t) jc?, if also Eqs. (3) and (18) are used. 

Eqs. (25)-(27) yield for the four-acceleration of the IDP 



~~d^ 



S 



pr 



C 



S 



mc u cr; 

J=2 



- 

tot 



C 



(28) 



In the further treatment we will consider the case a'p^j = for j = 1, 2. As a consequence, the 
equation of motion will be of the form 



dpi" _ 1 , 
1^ ~ c 



a LJ 



a. 



1 - 



c 



(29) 



and the four-acceleration will be 



dUi" 



Cr'pr S a CO 

mc u 



(30) 



The change of particle's mass is given by Eq. (27). Conventional approach is that the force due to the 
solar wind bombardment considers fixed mass of the particle (see, e. g., Mukai and Yamamoto 1982). 
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2.4. Total scattering cross section 

We have used the total scattering and pressure cross sections, in the previous theoretical parts. We 
need to determine them. We will consider spherical particle. 

We will use a some sort of approximation to the hard-core scattering problem, which corresponds 
to the limiting case of a short-range potential V{r) = oo for r < i?, V{r) = for r > R (hard-core 
potential, see, e. g., Iro 2002, p. 158). We will consider the scattering of point solar wind corpuscles 
from an almost hard sphere of radius R. In the case of the infinitely hard sphere of radius R, "the 
dynamics reduces to the laws of reflection at the surface of the sphere" (Iro 2002, p. 158). The result of 
classical physics is following: "In the case of a finite-range potential, the total [scattering] cross section 
is finite and gives the effective area of the potential. (This is actually the definition of a finite-range 
potential.) For example, when point masses are incident onto a hard sphere, ctj^j is the cross section 
of the sphere - only particles incident within that area are deflected." (Iro 2002, p. 161). 

However, correct physics for the incident electromagnetic radiation suggests that geometric cross 
section may not lead to correct results for the incident solar wind corpuscles (Klacka 2008a, 2008b). 
Inspiring by de Broglie's idea about wave character of massive particles, one; can come to the conclusion 
that scattering by a hard sphere at very high energies leads to the total scattering cross section 

ai,, = 2nR^ (31) 

and "the classical total cross section is just half of the quantum- mechanical result in the limit of very 

short wavelength" (Messiah 1999, pp. 393-395). 

If we use a some sort of approximation to the hard sphere, we can use the total scattering cross 
section given by Eq. (31), in our paper. 

As for the comparison of the results obtained by quantum and non-quantum physics, we will use 
Eq. (25) or Eq. (29). The non-quantum approach uses a^^j = a'p^. = A! = ttR"^: 

/ non— quantum \ / 

Comparison between Eqs. (29) and (32) yields (compare coeflicients at and J/'*): 

rr' — A' 

x' {quantum) = [x' {non — quantum) — 1] — ^ -|- 1 . (33) 



(7 



tot 



Using also Eq. (31), one obtains 

x' {quantum.) = ^ [x' {non — quantum) + 1] . (34) 

The case {o'pj. — A' , ct^qj — 2A'} is analogous to the cases of perfectly absorbing or refiecting 
spheres within geometrical optics approximation for electromagnetic radiation (Klacka 2008b). This 
analogy explains also the importance of quantum physics in our derivations - non-quantum physics 
would not yield correct results in the limit u — > c. 



2.5. Equation of motion to tiie second order in v/u 

In the approximation to the first order in v/c, we can replace the spacelike part of the four-acceleration 
of the IDP by acceleration dv/dt, where t is time measured in the stationary reference frame (associated 
with the Sun). Further, using Eqs. (2), (11), (14), (15), (17) and (22), we can express the right-hand 
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side of Eq. (30) in the approximation to the second order in v/u. We get 

1,2 



m 



V ■ u 



u 



\ V" ^ V ■ U V 

+ o — " + 

2 u u 



A' = itB? = cjL 



pr 



(35) 



where u = u/u is the unit vector in direction of the solar wind. 

Let us introduce the cylindrical coordinate system associated with th(^ orbital plane of the IDP 
and determined by unit vectors cr (radial vector), er (transversal vector) and bn = brx bt (normal 
vector). We can write (Klacka 1994) 



u 



7t = sin s 



where 

•jR = cose 
and 



ut = ^ k X Br = {bt cos i — bn cos O sin i) , 



(36) 



(37) 



= \j {cos if' + (cos 9)^ (sini)^ 



(38) 



The quantity e is an angle between the radial direction and the real direction of the solar wind. The 
unit vector k corresponds to the vector of angular velocity of solar rotation. The inclination of the 
orbital plane of the IDP with respect to the solar equatorial plane is i. Finally, is a position angle of 
the IDP (an angle measured from the ascending node of the orbit of the IDP to its actual position). 

Inserting Eqs. (36) and (38) to Eq. (35), and using the decomposition of the velocity vector into 
its radial and transversal components, v = vrbr + vtbt, one obtains 



dv 
'dt 

Xr 



A'nmiv? 



m 



{XrBr + XtBt - XnBn 



vr V cos I Vt V 

+ 1R 1- 7t 

u u N u u 



IR 

1 
1 



(1 + 7^) 



VR 



IRIT 



COS I Vt 



IR 


Vr 


- 7t 


cost 


Vt 


u 


N 


u 


IR 


Vr 


- 7t 


cosi 


Vt 


u 


N 


u 



N u 

1 v^ 

' 2 ^ 

1 f_ 

' 2 ^ 



1 v^ 

+ 7-R o — 

cosi 

7t 



7t 



N 

COS © sin i 
N 



Vt 
u 



(39) 



The angle e is small, its value lies between 2° - 3° (Bruno et al. 2003). Thus, we can neglect terms 
proportional to 7p and 7t(w/u)^. Similarly, we put jr ^ 1. Then 



dv 



A'nmiu^ 



1 v^ 



Vr 
u 

cosi 



cost Vt 



eR 



V u ) N 

vr\ cos 9 sin i 



Vt 
u 



Bt 



ejv + 



Vr 



u J 



(40) 
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3. Equation of motion - electromagnetic radiation effect 

Up to now, we dealt with the action of solar wind on the motion of an IDP. The role of solar electro- 
magnetic radiation cannot be neglected in the motion of the IDP in the Solar System. Relativistically 
covariant equation of motion for an arbitrarily shaped dust particle under the action of parallel beam 
of photons (Klacka 2008a, 2008b): 



d7 



(41) 



where is four-vector of the particle of mass m, four-vector of the world-velocity of the particle is 
given by Eq. (21) and four-vectors bj, j = 1, 2, 3 are given as: 



1^ 



.(»)(i-^) 



d = 



j{v) (1 + 



c 

V ■ e' 



V ■ e 



+ 



1, 2, 3 



(42) 



where {e'j-. j = 1, 2, 3} is orthonormal vector basis in the proper reference frame of the particle and 
{ej: j = 1,2.3} is corresponding vector basis in the stationary frame; ei corresponds to the radial 
direction (i.e. the Sun - particle direction). Seimg is flux density of the electromagnetic radiation and 
(7' (j = 1, 2, 3) are spectrally averaged cross sections of radiation pressure 



r' 

■ 



/o°°J(A) c;,.^ (A) dX 
/~/(A) dX 



i = 1,2,3, 



(43) 



where /(A) is the flux of monochromatic radiation energy. If (7^^ 2 — C'pr 3 = 0, then Eq. (41) reduces 
to the Poynting-Robertson effect (Poynting 1903, Robertson 1937, Klacka 2008a, 2008b, Klacka et al. 
2009), since also components of the thermal emission force j {j = 1, 2, 3) arc equal zero, in that 
case (Mishchenko 2001, Mishchenko et al. 2002). It can be easily verified that Eq. (41) yields dm/dr = 
0, i.e., mass of the particle is conserved, under the action of electromagnetic radiation. 
To the first order in u/c, Eq. (41) yields 



dv 



m c 



3 



+ [(1 + v)^^ - c. 



(1 - v-e'j/c)e'j + v/c , j = l,2,3. 



(44) 



It is worth mentioning to stress that the values of radiation pressure cross sections C'^^^^, j = 1, 
2, 3, depend on particle's orientation with respect to the incident radiation - their values are time 
dependent, in general. General equation of motion, represented by Eq. (41) or Eq. (44) differs from the 
Poynting-Robertson effect. Eqs. (41)-(44) hold for arbitrarily shaped particles. Experimental evidence 
that nonspherical dust grains move in a different way than spherical particles was given by Krauss 
and Wurm (2004). 
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As it is conventionally used in Solar System studies, we will restrict ourselves to the Poynting- 
Robertson effect, as for the effect of solar electromagnetic radiation. Thus, instead of Eq. (44), we will 

use 

at m c l\ c / c J 

where a dimcnsionlcss efficiency factor of radiation pressure i is defined by relation Q'^^ ^ = 
C'p^,i/A'- The values of Q^^^ I can be calculated according to Mie (Mie 1908; see also van de Hulst 
1981, Bohren and Huffman i983). 



3.1. Electromagnetic radiation effect as a special case of the solar wind effect 

The transformations 

i) w — > c ei, 

ii) S' = S aoj/u S'^i^g = Seimg wl (Klacka 2008b), 

iii) ^pr,j C'pr.j (i = 1> 2, 3), and, 

iv) x' = 1 (Klacka 2008b) 

reduce Eq. (25) into Eq. (41) without thermal emission terms {F^j = for j = 1, 2, 3). This means 
that our theory for electromagnetic radiation effect without thermal emission is consistent with the 
theory for corpuscular radiation effect. 



4. Equation of motion - solar radiation and solar gravity 

Let us consider spherical body orbiting the Sun under the action of solar radiation, i.e. solar corpuscular 
(solar wind) and electromagnetic radiation. The effect of the solar electromagnetic radiation on the 
motion of spherical particle corresponds to the Poynting- Robertson effect (P-R effect). 



4.1. Solar electromagnetic radiation effect 

It is useful to introduce a /3-parameter defined as the ratio of radial component of the radiation force 
and the gravitational force between the Sun and the particle with zero velocity: 

_ Le A' Q;, 

4 TT C TO /i 

II = G {Mq + m) = GMq . (46) 

Lq = 3.842 X 10^*^ W (BahcaU 2002) is luminosity of the Sun, Q'^^ = Q'^^ -^, G is the gravitational 
constant and Mq is mass of the Sun. For homogeneous spherical particle we can write 



Q' 

/3 = 5.760 xlO^ ^ 3-, (47) 

i?[/xm] p[kg/m ' 



where R is radius of the particle and p is mass density of the particle. Conventionally it is assumed 
that (i = const: neither optical properties nor mass of the IDP change. We do not restrict ourselves 
to the validity of this assumption. 

Now, using Eq. (46), the relation Seimg = -^0/(47rr^) holds, where r is a heliocentric distance of 
the IDP. On the basis of the decomposition of the velocity vector v = vrCr + vtGt, we can rewrite 
Eq. (45) into the form 



dtjp_ji r2 



1 - 2 

c 



-)en- - , (48) 
This is the acceleration of the IDP under the action of the Poynting-Robertson effect. 
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4.2. Solar wind effect 

Let us replace the fraction behind the curly braces in Eq. (40) by new quantities 

A'nmiu^ _ Tj ^ u jjL 
m ~ Q'^ c r 



Pz^- (49) 



T] is conventionally a constant of the value: 0.22 (Whipple 1967, Dohnanyi 1978) or 0.30 (Gustafson 
1994, Abe 2009). The solar wind speed u values of 350 km/s (Dohnanyi 1978) or 400 km/s (Gustafson 
1994) are used in modelling of orbital evolution under the action of solar wind. 

Let us look on the numerical values of r] and u on the basis of the solar physics data. In order 
to calculate the value of 77, wc need the values of n, u and Ei, on the basis of Eq. (15). The average 
values near the orbit of the Earth (1 AU) are (Hundhausen 1997, p. 92): proton density rii = 6.6 
cm~^, electron density n2 = 7.1 cm~^, He^"*" density = 0.25 cm~^, flow speed u = 450 km s~^. Eq. 
(15) yields for the average value S (solar wind; 1 AU) = u Y^^=i ^i-^i i = ^ = 515.642 

kg s~^. Moreover, we will take into account that riiU (1 AU) = (n,) (u) (1 AU) (1—0.15 cos(/5)^, (p = 2 
TT [t — t{max)]/T, T = 11.1 years and t{max) is the time of the solar cycle maximum (Svalgaard 1977 
- chapter 13). This result, together with Eqs. (11) ( w = 1), (14), (17) and (30), yields S{solar wind) 
= Seimg {A! ) a'^) T] and 

T] = r]o (1 - S cos(ff , 

u = uq (1 — 6 cos (fi) , 

m = 0.38 , 

6 = 0.15 , 

uo = 450 km/s , 

o t- tretard " t {max) 

V = 27r , 

T = 11.1 years , (50) 

if we put apr = A' and 5e;„g(lA[/) = Lq / [An [r{= lAU)]^ ], Lq = 3.842 x lO^^ W. The value of 
T represents the average value of the solar cycle period (see, e. g., Foukal 2004, p. 366). The retarded 
time tretard IS of the order r/uo and it is only a better approximation to reality than the omission of 
this term. 

4.2.1. More exact solution 

In reality, one needs to know concentration n(r, t) and solar wind velocity u{r, t), when the dust grain is 

situated at a position of heliocentric distance r at the time t. More precise information can be obtained 
from continuity equation (radial component of the velocity u is approximated by the magnitude u) 

dn 1 d(r^nu) ^ 

if r / 0. Using the observational fact n = const u / (Svalgaard 1977 - chapter 13), one obtains 

du dv^ „ 

if r 7^ 0. Using the boundary condition 

f „ I t — I (mux)'] 1 

lim u(r, t)=uo<l - 6 cos 2tt \ , (53) 

1 — '0 I 1 I 
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the quasi-linear partial differential equation Eq.(52) can be solved. 

Eq. (52) is known as Burgers equation: 
du/dt + u du/dx = , 
see, e.g., Sevcovic (2008, pp. 40 - 42). 
If the boundary condition 
u{x = 0, i) = x{t) 

is given, then the solution of the Burgers equation is 

m(x, t) = xit — x/u[x, t)). 

The last nonlinear algebraic equation can be solved by an iteration method, e. g.: 

u = x{t~ x/uq) + linife^oo Vk, 

■Vk+i = x{t - x/[x{t - x/uo) + Vk]} - x{t - x/uo), 

vi = 0. 

On the basis of the known solution of the Burgers equation, Eqs. (52)-(53) yield 



u{r, t) =uo<l — S cos 



27r 



t — r/ {2u (r, t)) — t{max) 



(54) 



Comparison with Eq. (50) shows that the retarded time is tretard = r/{2u). 

The nonlinear algebraic Eq. (54) can be solved by the iteration method presented above, or by the 
following iteration: 



Ufe+i(r,t) = uo\l - (5 cos 



2n 



t — r/ (2ufc (r, t)) — t{max) 



(55) 



since the right-hand side of Eq. (54) is a contractive/contraction function for the case u"^ > i: r 5 
I T and heliosphere is characterized by condition r < 150 AU (approximately). We can put u\(r,t) 
= Wo- 



4.2.2. Summary 

Using defition by Eq. (49), we can summarize, on the basis of Eqs. (50) and (54): 



m 


_ -q u II 


u{r, t) 


= uo 1 1 — ^ cos 


r]{r, t) 


= m[u{r,t)/uof , 


Vo 


= 0.38 , 


S 


= 0.15 , 


Uq 


= 450 km/s , 


T 


= 11.1 years . 



t — r/ {2u{r,t)) — t{max) 
2tt = 



(56) 

One can use also n = n{r, t) = uq [u{r, t)/wo] (1 AU / r [AU])^. These results represent a more realistic 
model than the model conventionally used. It takes into account more observational facts (Svalgaard 
1977, Hundhausen 1997). 

A more simple accesses use Eq. (50) with tretard = f / {2uq) or tretard = 0. Shock waves do not 
exist in these simple cases. The shock waves are generated by solution of the Burgers equation and 
this is considered in Eq. (56). As was pointed out in the comment to Eq. (55), the shock waves are 
not realized in the Solar System. 
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We put i = in Eq. (40). Then, the acceleration of the IDP caused by the solar wind has the form 



'dtj 



sw 



Q'vr ^ 



It 



!1 - 2^^ 
c 



2 — ) efl - 

c / 



Vt 



Vt 

— eR 
c 



1 

2 uc 



c 

Vr V 

c u 



V c c / 



(57) 



4.3. Equation of motion 

Gravitational acceleration from the Sun is —{ji/r'^)eR. Neglecting the solar wind pressure term 
P{ri/Q' ){iJ,/r'^){u/c) in Eq. (57), we can write the final equation of motion of the IDP in the form 



dv _ _ M (1 - /?) 
dt ~ r2 



eR 




+ 1^ — eR + — - \ ; 58 
2 uc c u] 

Eqs. (48) and (57) were used. We introduced a new central acceleration — /u(l — /3)eii/r^, i.e. the 
gravitational acceleration from the Sun reduced by the solar electromagnetic radiation pressure. Other 
terms on the right-hand side of Eq. (58) constitute the nongravitational disturbing acceleration. Eq. 
(50) has to be taken into account. 

Moreover, Eq. (27) can be rewritten to the form describing decrease of particle's radius R: 

dR K \u-v\ 2 

— = 2 (l-(5cos<p) , 

at u 

u = Uq {1 — 6 COSip) , 
5 = 0.15 , 
uq = 450 km/s , 

„ t- tretard ~ t (max) 
¥' = 277 , 

T = 11.1 years , (59) 

where K is a, constant characterizing decrease of the radius of the particle; Eq. (50) was used, also 
(see also text between Eqs. 50 and 57). Eq. (37) has to be used, too: u = u ii. One can use Eq. (59) 
as an approximation to the process of erosion of the particle due to the solar wind corpuscles. Since /? 
is a function of R (also Qp^. is a function of R), Eqs. (50), (58)-(59) have to be solved simultaneously, 
together with the Mie's calculations yielding Q'^^ for a given R. 

The real flux density of solar wind energy and the approximation of a constant flux (ry = rjo) for 
the radial solar wind, can be described by the approximation 

r) = r]o= 0.38 . (60) 



Analytical approach to solution of Eqs. (50), (58)-(59) is presented in Sec. 5. Detail numerical 
solutions are given in Sec. 6. 
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5. Secular evolution of particle's orbital elements under the action of solar radiation - 
analytical approach 

We have obtained complete equation of motion in the preceding section. We want to obtain qualitative 
understanding of the orbital evolution of IDP. This task is the main subject of this section. Trend 
is that /3 is a decreasing function of R. Thus, the term — [/x(l — f})/r^] br, in Eq. (58), does not 
correspond to Keplerian acceleration, if Eq. (59) is taken into account. As a consequence, osculating 
orbital elements have to be calculated for Keplerian acceleration and it is given by the term — (/x/r^) 
eR. 

We will calculate the secular evolution of semimajor axis a, eccentricity e of the particle's orbit 
under this nongravitational perturbation. We have to use Eqs. (100) or (103) in Klacka (2004) and 
Eqs. (32), (34) and (37) in Klacka (1993b), for the time evolution of the osculating elements. 

Secular values of semimajor axis and eccentricity, when central acceleration is given by solar gravity 
term — (/x/r^) br, are 



1 



/? ^1 + e| + 2e/3 cos x 



(1 + e/jcosx) 



dx 



1 /-a^ J(l-/?)^e^+/32-2/3(l-/J)e^cosa; 



(1 + 6/3 cos x) 



dx 



(61) 



The quantities denoted by index /? hold for central acceleration — [;u(l — P)/r^] br and it is supposed 
that their values do not significantly change during particle's revolution around the Sun. Analytical 
approach for secular evolution of orbital elements is possible under such an assumption, otherwise 
detail numerical integration of Eq. (58) has to be done. The quantities denoted by index /? can be 
calculated from Eq. (67) below. 



5.1. Calculation of quantities for Eq. (61) 

Let us use perturbation equations of celestial mechanics in the following form: 

a/3 



da/} 
~dF 



dep 
dt 



dt 




+ 



/i(l-/3) 
/3 



ttR eissmffj + aT(l + C/jcos//?)] 



e/3 + cos fp 



a,R sin //3 + ar cos + 



1 + 6/3 cos fp 



/3 



{ep + cos fp) , 



1 

+ — 



1 

(3) ep 



ClR COS fp 



. 2 + 6/3 COS fp 

arsmfp - ^ 



ep cos fp 



ep l-p 



sin//3 , 



(62) 



where pp = ap{l — e'p), fp is true anomaly of the IDP and ujp is argument of perihehon of the particle's 
orbit; the dot over $ denotes differentiation with respect to time. It is assumed that the longitude of 
the ascending node is time independent. Moreover, ur and ar are radial and transversal components 
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of the disturbing acceleration. Eqs. (62) are consistent with Klacka (1993a - Eqs. 12, 14) and Klacka 
(1993b - Eqs. 32, 34, 37 if M = Mq{1 - 13)). Using Eq. (58) and expressions 



vr = 
we obtain 



hii-P) 



P0 



e/3 sin ffs , vt = 



Ifi{l-f3) 



P0 



(1 +6/3 cos /^) , 



(63) 



an 



2 1 



V \ o 1^ (1 - P) IP0 . r 

Q^r^ — -c — ^^^^^^'^ 



m(1-/3)/P/3 



7T 



P/3 



(l + e/3 cosfp) 



+ 



\ {I + el + 2ef} cos ff}) + siv? ff} | 



At 



1 + 



Q 



P (1 + e/3COS//3) 



pr 

J.2 f. 



+ 



Q'pr 

m(1-/3)/P/3 



7t u 



It 



6/3 sin ff} 



6/3 (1 + 6/3 COS //s) sin//3 



(64) 



Secular evolution of the orbital element g we get by the time averaging of dg/dt over one orbital period 
P, i.e. 



dg 



dt 



1 



^ = ^ -^dt 



dg 



dt 



(65) 



We have used the second and the third Kepler's laws: r^dfp/dt = — f3)pi3 — dtop/dt — {dftp/dt) 

cosi/3 = ^ /x(l — P)p/3 and a^/P'^ = /u(l — /3)/(47r^); p/3 = a/3(l — e^), a;/3 is argument of perihelion 
and f2/3 is longitude of the ascending node. 

As for the quantity $, it is a function of radius R of the particle. Using an approximation (3 = 
A/R + B, where A and B are constants for a given particle, we can write 



dp _ A dR _ A K 
m ^ ~ R^ ~dt ^ R^ ^ 



(66) 



if also dominant part of Eq. (59) is used {6 = 0, |ti — t;| = u). 

Application of Eq. (65) to Eqs. (62) and (66), using the assumption that orbital elements and 
/3-paramctcr do not significantly change during particle's revolution around the Sun, yields (5 = is 
assumed) 



dap ^ ^ 1^ 
dt c 



2 + 3el 



0-/3 



{l-ei) 



3/2 



Qpr 



27T 



Uo 



2 + 3e2 ^f,{l-p)/p^ 
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dt 



dt 



1 + 4 



1-13 B?^^ (. 2\3/2 ' 



II 5e/3/2 



^ awl e| 



X M + 



0; 



1 - gTT 



e/3 



1 - p2 



1 



J2_ /3 ^ . 



X < 7T 



73 



1 - p2 



1 

2 uo 



(67) 



where the symbol ( ) was omitted on the left-hand sides. The quantities with index (3 in Eq. (67) are 
the quantities which have to be used in Eq. (61). Eqs. (67) hold under the assumption r] = rjo. 



5.2. Complete set of differential equations for secular evolution of particle's orbital elements 

Wc have obtained the following set of differential equations for secular evolution of orbital elements: 
Eqs. (61), (66)-(67). This set has to be completed by equation for secular evolution of radius of the 
particle. Using Eqs. (59) and (65), we obtain {5 = 0, \u — v\ = u) 



dR 
'dt 



K 



al^l 



(68) 



Initial conditions must be added to the set of differential equations for orbital elements. If the 
particle is ejected from a parent body of known orbital elements, then the particle's initial orbital 
elements have to be calculated from Eqs. (60)-(61) in Klacka (2004). 

If one would take into account also thermal change of optical properties of the spherical dust 
particle, then also further changes of orbital elements exist (secular changes of semimajor axis and 
eccentricity, perihelion motion). This would correspond to the change of parameter A and B in Eq. 
(66). We do not deal with this case (see Klacka et al. 2007, Pastor et al. 2009). 



5.3. Discussion 

Let us look on secular evolution of semimajor axis of the particle's orbit. If the P-R effect and 
radial velocity component of the solar wind are considered alone (i.e. 77- = and /3 = 0), then Eqs. 
(61) and (66)-(68) show that a is a decreasing function of time. If one takes into account also the 
non-radial velocity component of the solar wind, or, /? > 0, the situation may be different. The secular 
value of a can be an increasing function of time. Thus, the effect of real solar wind may cause particle's 
spiralling outward from the Sun. (Similarly, also the secular value of e can be an increasing function 
of time.) 
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5.4. Radial solar wind and decrease of particle's radius 

According to Eqs. (27) and (59), the mass of the particle may decrease. Eq. (68) holds under assump- 
tion that the decrease of the particle's radius is small enough during particle's revolution around the 
Sun. 

When we consider only the P-R effect and radial solar wind effect (7t = 0), we are able to calculate 
the radius of IDP as a function of its orbital eccentricity. Using Eqs. (67) and (68), one can immediately 
write 




dR ^ [2 \^ Q'p^J c 1-/3 i?2 
Using also Eqs. (47) and (66), Eq. (69) can be integrated: 

6/3 = epin < —. ^r-^ ^ ^xp [ki {R - Ri„)] 



(69) 



^ A-{l-B)R \Rin 
5 iJ, B 
2 c K[cm AU^ yr-i] 

5 m 1 J . 5.760 X 102 



'^=2ci^[cm AUVr-V^^''' 
IJ, = Att^AU^ yr-2 , 



kg/m^ 



c 



6.3114 X 10'' AU yr-i , 



V = 0.38 . (70) 

and the quantities A and R are given in microns; the subscript "in" denotes initial values. 

In order to find secular evolution of semimajor axis a and eccentricity e, according to Eq. (61), we 
have to solve the following set of equations: Eq. (66) for a given values of A and B, equation for 
in Eq. (67), Eq. (68) and Eq. (70) (or, equation for C/j in Eq. 67 instead of Eq. 70). 

5.5. Semi-latus rectum, time of spiralling for radial solar wind and K = 

Prom Eqs. (67) we can determine also the secular evolution of the semi-latus rectum of the particle's 
orbit. Whence = ap{l — e|) we can write 

what using Eqs. (67) yields (the symbol ( ) is omitted) 

/ \ 3/2 

dpfj 
dt 




c 1-/3 i?2 




C ^fi{l-p)/p0 pp 



1-4 

^ (72) 
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Let us rewrite equation for {dcfj/dt) into the form (the symbol ( ) is omitted) 

3/2 



dt 



e/3 



('- 



4)' 



Q vr 



A 



l-(3 B? 



K 



Uo 



Q'pr C ^^(l-/3)/p0 



4) 



3/2 



P} 6/3 



(73) 



Now, let us consider only the P-R effect and the radial solar wind effect. We put 7^ = in Eqs. 
(72)-(73). In this case, we obtain the following equation 

de/3 6/3 ^ ^ ^ 

f / \ ^ -1 

I 1^ / \ ,, 1 \ \ 

(74) 




from Eqs. (72) and (73). Eq. (74) yields the relation 

4/5 



V0 = Pf} 



6/3 

6/3 ir> 



K = Q 



(75) 



where pfj in and e/3 m are initial values of semi-latus rectum and eccentricity of the particle's orbit. 
Eq. (75) can be considered as a generalization of the result obtained by Wyatt and Whipple (1950): 
we have taken into account not only the P-R effect, but also the radial solar wind effect. Eq. (75) 
allows us to write equation for secular evolution of eccentricity in the form 



\dt / 



1 + 



Q'm 



8/5 
M ^13 in 



M) 

3/5 



3/2 



K 



. 



(76) 



It is evident from Eqs. (74) that due to the P-R effect and radial solar wind the particle is spiralling 
inward to the Sun, for K = 0. Semimajor axis and eccentricity e/3 of the particle's orbit converge 
to (see Eqs. 75, 76 and relation among a/3, Pjs and e/3). Eq. (76) can offer the time of spiralling of 
the particle with initial orbital elements a/3 ,„ and e/3 i„ into the orbit with osculating elements a/}, 
ep. This time is given by relation 



T(e/3 in, 6/3) 



1 + 



Q' 



^0 in 



(1 ~ 6^ m) 



pr . 



e/3) 

K 



n X 

Je, (1 - 



8/5 
^/3 in 



I {ep in,ei3) , 



,3/5 



'e0 

. 



2^3/2 



dx 



(77) 
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The time of spiralling of the particle into the Sun is then 



T (e/3 m,0) = - g 



2 

r / \ -1-1 / 1 ^2 \ 

2 




in 8/5 m)0) ) 

in 

X = . (78) 



Let us consider two particles characterized by the values (3i and ^2- Moreover, let the particles 
have the same value of Q'pj.- If we are interested in times of stay of the particles within an interval of 
semimajor axes {aiower, dupper), then Eq. (77) yields: ti/t2 = (32/ if the initial values of semimajor 
axes and eccentricities are equal for both particles. On the basis of Eqs. (72)-(73) {K = 0), this result 
can be approximately generalized also to the case of real solar wind effect, under the assumption 
/32 « 1. 

As for the secular evolution, under the assumption that particle's radius does not decrease, we 
have to solve only one differential equation Eq. (76); Eq. (75) immediately yields semi-latus rectum 
and one can easily obtain semimajor axis ap = / {1 — e|). 

5.6. Summary 

Our analytical approach shows that decrease of particle radius due to the solar wind abrasion (cor- 
puscular sputtering) generates increase of particle's semimajor axis and eccentricity. This result is 
consistent with detailed numerical calculations presented by Kocifaj and Klacka (2008). 



6. Numerical results 

In this section we will concentrate on orbital evolution of interplanetary dust particle under the action 
of solar electromagnetic and corpuscular radiation when solar wind erosion is neglected. The results 
based on our new approach presented in Sees. 2-4 will be compared with the standard approach when 
only radial solar wind with constant r] is taken into account. 



6.1. Radial solar wind 

Real flux density of solar wind energy and the approximation of a constant flux (77= 1/3) are compared, 
for radial solar wind, in Fig. 1. Eqs. (50) and (58) were numerically calculated for a particle of 
fixed radius: /3 = 0.01, Q'^^ = 1. The orbital elements were calculated for the central acceleration — 

(GMe/r^) br. 

The variable solar wind properties cause that semimajor axis exhibits a little faster decrease than 

the standard approach when 77 = 1/3 is used. The same holds also for eccentricity of the particle. The 
time of spiralling toward the Sun for variable rj is about 10 % smaller than for the case of = 1 /3. The 
most significant difference between variable and constant // exists for secular evolution of argument 
of perihelion. While the constant rj produces no shift of perihelion, the more realistic approach of 
variable rj produces nonnegligible shift of perihelion. Moreover, the shift of perihelion depends on the 
initial position of the particle in its orbit. Fig. 1 shows evolution of the shifts of perihelia for two 
cases, for aphelion and perihelion ejections of the particle from a parent body (zero ejection velocity). 
Semimajor axis and eccentricity of the parent body were ap = 5 AU, ep = 0.9. 



6.2. Real solar wind 

If we take into account even more realistic description of the solar wind, when its nonradial velocity 
component is considered, the resulting orbital evolution differs from the cases discussed in Sec. 6.1. 
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t [10* years] 

Fig. 1. Orbital evolution of spherical dust particle (Qp^. = 1, /? = 0.01) for radial solar wind. Results 
for variable rj and standard value ry = 1/3 are compared. Secular evolution of semimajor axis and 
eccentricity and also shift of perihelion are depicted. More realistic case yields that the evolution of 
the shift of perihelion depends on the position (A - aphelion, P - perihelion) of the parent body at the 
time of the particle's ejection. Effect of solar wind erosion is of negligible importance. 

Fig. 2 compares the standard approach (radial solar wind, ry = 1/3) with the variable r] and jt 0. 
Again, Eqs. (50) and (58) were numerically calculated for the particle of fixed radius: /3 = 0.01, Qp^. 
= 1. The orbital elements were calculated for the central acceleration — (GMo/r^) en. 

The real solar wind causes that semimajor axis exhibits a little slower decrease than the standard 
approach (radial wind, ij = 1/3) is used. The same holds also for eccentricity of the particle. The time 
of spiralling toward the Sun for real solar wind action is about 15 % greater than for the case of 77 = 
1/3 (one must be aware that this value holds for the case presented in Fig. 2 - the real percentage may 
be greater/less for larger/smaller initial values of semimajor axis than for the initial 5 AU). The most 
significant difference between the usage of the real solar wind and the standard approach exists for 
secular evolution of argument of perihelion. While the constant r] produces no shift of perihelion, the 
realistic approach produces nonnegligible shift of perihelion. Moreover, the shift of perihelion depends 
on the initial position of the particle in its orbit. Fig. 2 shows evolution of the shifts of perihelia for two 
cases, the aphelion and perihelion ejections of the particle from a parent body (zero ejection velocity). 
Semimajor axis and eccentricity of the parent body were ap = 5 AU, ep = 0.9. 

7. Discussion 

We have derived relativistically covariant equation of motion for the action of solar wind corpuscles on 
motion of interplanetary dust particles. As for spherical shape of the particles, the equation of motion 
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radial SW - 











2 4 6 

t [lo" years] 



Fig. 2. Orbital evolution of spherical dust particle, Q'p^. = 1, /? = 0.01. Secular evolution of semimajor 

axis and eccentricity and also shift of perihelion are depicted. Results for real solar wind (variable 77 
and non-radial component of solar wind velocity vector s = 2.9° are included) and radial solar wind 
(standard value r] = 1/3) are compared. More realistic case yields that the evolution of the shift of 
perihelion depends on the position (A - aphelion, P - perihelion) of the parent body at the time of 
the particle's ejection. Effect of solar wind erosion is of negligible importance. 



is represented by Eq. (29). It differs from the force conventionally presented in literature (although 
only to the first order in v/u): 

Fs^ = [ (1 - 2 r/u) ER - {re/u) bt ] , 

where v is the velocity of the grain v = r en + r0 et, u is the heliocentric solar- wind speed and F^^ 
is the force on the dust for u = (Minato et al. 2004 - Sec. 2.1; equivalent force is presented in Eq. 
(7.10) by Mann 2009; see also p. 12 in Burns et al. 1979). 

We have to stress that the standard form corresponds to Eq. (32) if a;' = 1 (reality is: 1 < a;' < 
3, approximately) and not to physically correct form given by Eq. (29). Eq. (29) yields the correct 
limiting result u ^ c equivalent to the P-R effect. 

As for practical application of the above discussed physical results, Eqs. (56)-(59) are astronom- 
ically relevant. Decrease of mass of the particle may cause its spiralling outward from the Sun and 
not toward the Sun, as it is commonly accepted for the action of solar wind on interplanetary dust 
particles. This result is evident from analytical equations presented in Sec. 5 (see, e. g., Eqs. 61-62) 
and detailed numerical calculations confirming this result can be found in Kocifaj a Klacka (2008). 
Sec. 6 concentrates on the action of solar wind on motion of interplanetary dust particles for the 
cases when variable flux of solar wind energy and non-radial solar wind velocity are considered. The 
most evident difference between the standard and more realistic approaches is represented by the shift 
of perihelion (secular evolution of argument of perihelion). It is generally believed that the shift of 
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perihelion does not exist, both for the solar wind effect and the P-R effect. However, the real action 
of the solar wind differs from the action of the P-R effect - the P-R effect really produces no shift 
of perihelion. The real shift of perihelion depends on the initial orbital position of the particle. The 
non-radial solar wind velocity can lead to outspiralling from the Sun, in the region of outer planets. 
The particle may or may not spiral toward the Sun due to the simultaneous action of the P-R effect 
and solar wind effect (see also Klacka et al. 2008). 



8. Conclusion 

Relativistically covariant equation of motion for arbitrarily shaped dust particle under the action 
of solar wind is derived. Change of the particle's mass is an indispensable part of the space-time 
formulation of the equation of motion for the action of the solar wind. The solar wind effect would 
reduce to the Poynting-Robertson effect, in the limiting case when: i) solar wind speed would tend to 
the speed of light, ii) no decrease of mass of the interplanetary dust particle would exist, and, iii) the 
velocity of the solar wind would be radial. However, the solar wind may have qualitatively different 
effect on orbital evolution of interplanetary dust particle, since the points ii) and iii) are not fulfilled, 
in general. The decrease of mass of the interplanetary dust particle and non-radial component of the 
solar wind velocity may cause outspiralling of the particle from the Sun. Time variable solar wind 
leads to the shift of perihelion of the particle. The found results may have important consequences for 
evolution of dust disks in the vicinity of stars with stellar winds. 



Appendix A: Emission from the particle 



The other possible force influencing dynamics of dust particle may originate from an emission, e.g., 
radioactive decay. Let the particle emits an energy E'^^ per unit tinic^ due to the emission in the 
proper reference frame. We will suppose that this emission is represented by the flux of particles with 
a speed u'^^. Further, we declare the orthonormal vector basis {f'j; j = 1, 2, 3}, as it was used in 



Sec. 2.2. The corresponding velocities are: u'e 



= u' 



J'j, j = 1, 2, 3. 



The outgoing four-momentum of the emission per unit time, in the proper reference frame of the 
particle, is 



\\ E' ■ - E' V r' ^^LJ^Uhl] 
1 c ^™ ' c f^^ ' c I 



(A.1) 



where r'j {j = 1, 2, 3) are dimensionless coefficients expressing the part of the total flux of radiation 
which is emitted in the corresponding directions. 

Lorentz transformations of (A.l) yield the outgoing four-momentum per unit time in stationary 
reference frame: 



Ie' ^ + ^E' r' - ^] 

r em ^ ~ f, eTO /-^ ' j y^emj ^ J 



(A.2) 



where 



= 7(«)(1 , 



Ue 



J = 



1, 2, 3 . 



V ■ u' 



u 



(7M-1) 



V ■ u 



(A.3) 
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Thus, expression on the right-hand side of Eq. (A. 2) should be added to the right-hand side of Eq. 
(21) when we would like to take into account also the effect of the emission of dust particle. 
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